Growth of a vortex polycrystal in type II superconductors 
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We discuss the formation of a vortex polycrystal in type II superconductors from the competition 
between pinning and elastic forces. We compute the elastic energy of a deformed grain boundary, 
that is strongly non-local, and obtain the depinning stress for weak and strong pinning. Our es- 
timates for the grain size dependence on the magnetic field strength are in good agreement with 
previous experiments on NbMo. Finally, we discuss the effect of thermal noise on grain growth. 



Understanding the phase diagram of high temperature 
superconductors is still a formidable challenge of mod- 
ern condensed matter physics. Typically high T c ma- 
terials behave in a magnetic field as type II supercon- 
ductors, with further complications due to the broader 
phase space — in terms of temperature T and field H — 
in comparison to conventional superconductors 0, 0, • 
Raising the temperature, the Abrikosov vortex lattice |j| 
melts into a liquid 0,1a, Li, while quenched disorder leads 
to more complex phases such as the vortex glass @ or 
the Bose glass Q. The natural question posed to the 
theorist is to explain the occurence of the various vor- 
tex phases, linking the experimentally observed behavior 
with the material microstructure. 

The Bitter decoration technique provides a powerful 
method to investigate the geometrical and topological 
properties of vortex matter by direct imaging of the vor- 
tices [10| . Its application to conventional superconduc- 
tors provided the first direct evidence of the vortex lattice 
fl"ll | predicted by Abrikosov 0. The lattice structure is 
often observed to coexist with topological defects, such 
as isolated dislocations, dislocation dipoles and grain 
boundaries. These last extended defects are the signature 
of a vortex polycrystal with crystalline grains of different 
orientations }wL Il2 ] . Vortex polycrystals have been ob- 
served, after field cooling, in various superconducting ma- 
terials such as NbMo [13111, NbS 2 pHHH^, BSSCO 
[3 and YBCO [HI. The grain size is typically found to 
grow with applied magnetic field [r3 . [T^ |. Moreover, two- 
sided decoration experiments show that the grai n bound- 
aries thread the sample from top to bottom |13L . i.e., 
one observes a columnar grain structure. Despite the 
wealth of experimental observations, there is no detailed 
theory accounting for the formation of vortex polycrys- 
tals. The issue is particularly interesting since recent 
experiments indicate that the reentrant disordered vor- 
tex phase of NbS2, comm only believed to be amorphous, 
is instead polycrystalline Whether this is a stable 

thermodynamic phase is still an open question. 

The theoretical description of vortex matter in high T c 
superconductors is centered on the role of quenched dis- 
order. Early theoretical considerations seemed to imply 



that even a small amount of disorder would lead to the 
loss of long-range order and to the formation of an 
amorphous vortex glass phase Experimental obser- 
vations did not confirm this view, since ordered vortex 
structures are typically observed even in presence of dis- 
order. This contradiction was resolved by a more detailed 
theoretical analysis of the weak disorder limit, show- 
ing a topologically ordered, long-range correlated phase, 
termed the Bragg glass [l9|. While the presence of a 
Bragg glass has been confirmed experimentally [2(j, the 
precise nature of the transitions into the amorphous and 
liquid phases is still debated. Recent theories highlight 
the importance of dislocations as mediators of the transi- 
tion [21| , in contrast with more traditional melting theo- 
ries, based on the Lindeman criterium lj . The properties 
of dislocations in the vortex lattice have been the object 
of extensive theoretical investigations p^ . l23ll23.l25| . but 
less is known about grain boundaries. 

Here we address the problem of the formation of a vor- 
tex polycrystal from the point of view of grain growth. 
In field cooling experiments, magnetic flux is already 
present in the sample as it is quenched in the mixed 
phase. Thus it is reasonable to expect that vortices are 
originally disordered and that, due to their mutual in- 
teractions, undergo a local ordering process through the 
growth of grains with various orientations which, in turn, 
implies the annihilation of several dislocation lines as well 
as their organization into grain boundaries separating the 
crystalline grains (see Fig. ^ for an example). The effect 
of quenched disorder is to pin the grain boundaries, hin- 
dering the growth process. Thus to understand the prop- 
erties of vortex polycrystals, we analyze the dynamics of 
grain boundaries in vortex matter as they interact with 
disorder. 

We first use continuum anisotropic elasticity to evalu- 
ate the elastic response of an extended three dimensional 
grain boundary to small perturbations. In the large- 
wavelength limit, the grain boundary elastic energy is 
found to be strongly non-local, with a linear wave vector 
dependence in Fourier space. Hence, it is not feasible to 
use a surface-energy approximation similar to the line- 
tension approximation frequently used to describe the 
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FIG. 1: Grain boundary structure from a simulation of inte- 
acting vortices after a sudden field cooling from a disordered 
vortex state. The colored non six-fold coordinated vortices 
indicate dislocations, mostly arranged into grain boundaries. 



elastic response of isolated dislocations. We thus use the 
non-local elastic energy in the framework of weak and 
strong pinning theories and estimate the grain size of 
a vortex polycrystal. Our results are in good agreement 
with the experiments, reproducing the dependence of the 
grain size on the applied magnetic field found in NbMo 
single crystals. In addition, we discuss the effect of ther- 
mal activation and the associated creep laws. 

The elastic energy of the vortex lattice can be ex- 
pressed in terms of the vortex displacement field u 

H=~ y f i 3 r[ Cll (Vu) 2 + (c 11 -c 66 )(V-u) 2 + C44(c),u) 2 ], 

where en, C44 and cqq are respectively the compres- 
sion, tilt and shear moduli, and we assume the applied 
field to point along the z direction. In the simplest 
continuum and non-dispersive approximation the elas- 
tic moduli can be estimated as en ~ C44 ~ _B 2 /47r and 
C66 — $oB/(8ttA) 2 , where B is the magnetic induction, 
$0 is the magnetic flux quantum and A is the London 
penetration length. 

We describe an extended grain boundary as an infinite 
array of edge dislocations arranged regularly along the y 
axis, with spacing D and Burgers vector in the x direc- 
tion. Note that due to the columnar grain structure of 
the vortex polycrystal, tilt and mixed boundaries which 
would require a more complicated description in terms of 
screw and edge dislocations, cannot occur. To obtain the 
elastic response of the grain boundary to small perturba- 
tions, we consider a generic deformation v = v n (z)x for 
the nth dislocation along its glide plane and find the elas- 
tic displacements solving the elastic equation associated 
with Eq. (JTJ with the appropriate constraints induced by 
the dislocations. The grain boundary elastic energy can 
then be obtained as in Ref. [24[ by a suitable expansion 
of Eq. QJ. The full calculation is somewhat involved and 



will be reported in detail elsewhere. Here we just quote 
the resulting elastic energy which in Fourier space reads 
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kernel is given by 



M(k y ,k z ) = 2c 66 \k y \ 
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The long-distance behavior of the kernel is captured by 
the behavior at small wave vectors that is given by 



M(k y ,k z ) ~ 2c 66 \k y \ + y/c 44 c 66 \k z 



(4) 



corresponding to a long-range non-local interaction in 
real space. It is convenient to work in an istropic refer- 
ence frame, rescaling the y coordinate by a factor \ 

The elastic kernel then becomes M(k) ~ K\k\, 
with K = y / C44C(36. In thin films, we can neglect the 
deformations along z and the kernel is simply given by 
M{k) ~ 2c66|fc|. Note that for an isolated vortex lattice 
dislocation it is correct up to a logarithmic factor to ap- 
proximate the elastic energy by an effective line tension 
[24j, but a similar procedure is not possible for a grain 
boundary. 

Quenched disorder induces elastic deformations of the 
grain boundaries and the competition between elasticity, 
disorder and a driving force acting on the boundary can 
be analyzed in the framework of pinning theories 0, [l8| . 
Driving forces for grain boundary motion can be exter- 
nally induced by a current flowing in the superconductor 
|25j| or, as in our case, internally generated by the or- 
dering process during grain growth [2||. Grain growth 
is driven by a reduction in energy: For an average grain 
size R and straight grain boundaries, the characteristic 
energy stored per unit volume in the form of grain bound- 
ary dislocations is of the order of Tq / R, where Tq is the 
energy per unit area of a grain boundary. Hence, the 
energy gain achieved by increasing the grain size by dR 
is T /R 2 dR. 

Physically, the removal of grain boundary dislocations 
occurs through the motion of junction points in the grain 
boundary network. As junction points must drag the 
connecting boundary with them, which may be pinned 
by disorder, motion can only occur if the energy gain at 
least matches the dissipative work which has to be done 
against the pinning forces. The dissipative work per unit 
volume expended in moving all grain boundaries by dR 
is a c b/ (DR)dR, where a c is the pinning force per unit 
area. Balancing against the energy gain yields the limit 
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In order to obtain an explicit expression for the grain size 
we distinguish between weak and strong pinning regimes. 

Weak pinning: In this regime, the grain boundary in- 
teracts with random stresses a xy induced by the elas- 
tic deformations of the vortex lattice |2lL l25j . The ran- 
dom stress correlations can be obtained directly from the 
vortex displacement correlations in the random manifold 
(RM) and Bragg glass (BG) regimes [21| 
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where R a is the length at which the vortex displacements 
become of the order of the lattice spacing a and corre- 
sponds to the onset of the BG regime. The roughness 
exponent in the RM is estimated as £rm ~ 1/5. The 
pinning energy for the grain boundary is then given by 
Hpin — J2 n I dzv(nD,z)b<j xy [v,nD,z\. The roughness 
of the grain boundary can be obtained by a scaling ar- 
gument [25| comparing this expression with the elastic 
energy in Eq. (0) , yielding Cgb — Crm in the RM regime 
and a logarithmic roughness for R > R a . Notice that 
both isolated dislocations and "dislocation bundles" are 
found to be rougher than grain boundaries [25|. which 
implies that the latter are more stable. 

The depinning stress can be computed within the 
framework of collective pinning theory: the energy asso- 
ciated with bending a grain boundary fraction of linear 
dimension L < R a over the characteristic distance v can 
be estimated as 
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where the first term represents the elastic energy, the sec- 
ond the pinning energy and the third the work done by 
an external driving stress a in displacing the boundary. 
Minimizing the first two terms, for v — a ~b |25| . we ob- 
tain the "plastic" Larkin length L p ~ (b/D) 5 R a , which is 
typically smaller than R a . The depinning stress is identi- 
fied as the stress necessary to depin a section of dimension 
L p : a c — Kb 2 /(DLp). Combining this expression with 
Eq. ©, using T ~ Kb 2 /D, we obtain R g — R a . The 
identification of R g with R a was proposed in Ref. [l2^ . 
but was not confirmed by experiments (see Fig. [21 and 
Ref. 01) E3- We therefore propose to interpret the ex- 
perimental data under a strong pinning assumption. 

Strong pinning: In this regime, pinning centers are 
strong and localized and one can assume that the dis- 
locations forming the grain boundary are pinned by in- 
dividual obstacles. We consider here the case of colum- 
nar defects, oriented along the z axis. This case should 
be relevant for the experiments of Ref. |l2j where grain 
boundary pinning is provided by screw dislocations in 



the superconducting crystal. In these conditions, the 
problem becomes effectively two dimensional (2D) and 
we can directly generalize the strong pinning theory of 
Friedel (2^, developed in the context of single disloca- 
tions, by taking into account non-local elastic interac- 
tions. The basic idea is to consider a grain boundary 
segment as it depins from a pair of strong obstacles. 
The length L of the segment corresponds to the effec- 
tive spacing between obstacles along the grain bound- 
ary, and it forms a bulge of maximum width v. After 
the grain boundary segment overcomes the pin it will 
travel by an amount which is, again, of the order of v 
and, hence, sweep an area of the order of Lv. At the 
depinning threshold, the grain boundary starts to move 
through a sequence of statistically equivalent configura- 
tions, and the freed segment will encounter, on average, 
precisely one new obstacle in the course of this process. 
This argument leads to the condition Lv ~ l/p where 
p is the area density of pinning defects. The elastic en- 
ergy per unit length of the bulge of width v and exten- 
sion L is 2c 66 6 2 i> 2 /D 2 (2D result), and should balance 
the work per unit length abLv/D done by the driving 
stress a in bowing the boundary. This energy balance 
provides a relation between L and v. Furthermore, at 
depinning the total force dabL/D should be equal to 
the defect strength /o, where d is the sample thickness. 
Combining the equations above we obtain the depinning 
stress <j c b = Dfo/(dLf), where the Friedel length Lf is 
given by Lf = 2c§§b 2 dj (fopD 2 ). Inserting the expression 
for the critical stress in Eq. JSJ together with the scale- 
independent surface tension Tq = 2c§§b 2 / D, we obtain 
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In order to use this result to fit the data in Ref. [T^| . 
we have to express it in terms of the reduced field 
B = B/H C 2, where H C2 is the upper critical field of 
the superconductor. The field dependence is implicit in 
the parameters b and D, i.e. b ~ D ~ a ~ B^ 1 / 2 , 
as well as in the shear modulus cqq ~ B, and in the 
pinning strength Jo- The pinning force due to a screw 
dislocation was computed in Ref. p9l | and is given by 

f oc BVa(i_B)ln(e/2.76oB) « B^ln(^/2.76 OJ B) o , 
where £ ~ 100A is the coherence length |3(|, and bp ~ 5A 
is the Burgers vector of the screw dislocation 29] . The 
resulting expression predicts a linear field dependence of 
the grain size with logarithmic corrections. In Fig. [21 
we can corroborate that the agreement of this prediction 
with magnetic decoration data from Ref. [l2T ] is quite sat- 
isfactory, especially if compared to the estimate based on 
local elasticity assumptions. 

In the discussion above we have neglected thermal fluc- 
tuations, which could induce an activated motion of the 
grain boundaries, particularly in high T c materials. This 
problem can be approached generalizing scaling theories 
of creep for vortices and dislocations 0, |25| , In the 
weak pinning regime, the relevant energy barrier that 
the grain boundaries have to surmount under an applied 
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FIG. 2: The grain size of a vortex polycrystal experimen- 
tally obtained from Bitter decoration of a NbMo sample as a 
function of the applied magnetic field (Ref. ^3) i s compared 
with the theoretical predictions. The calculation based on 
Friedel statistics with non-local elasticity compares favorably 
with the data. For comparison we report as well the result 
obtained using Friedel theory with local elasticity and the ear- 
lier estimate from Ref. |T3| . formally equivalent to our weak 
pinning result. 



stress a < a c is given by U(a) = f/o(o' c /cr) A ', where 
Uq ~ Kb 3 R a and fj, = 1. In our case, the applied stress 
is the ordering stress, so that we have <r c /a ~ R/R a . 



Using this expression in the energy barrier for thermally 
activated grain growth, it follows 



dR 

t— = i? a exp 
at 



kT R a 



(9) 



where r is the appropriate characteristic time. The equa- 
tion can readily be solved yielding, in the long time limit, 
a logarithmic growth R(t)/R a = fcT/£/n log(t/V). This 
law holds for R > R a when the grain boundaries would 
be pinned at T = 0. In the initial growth stage R< R„, 
we can neglect pinning forces and the dynamics is ruled 
by the ordering stress: R ~ 1/-R, yielding a power law 
growth R(t) ~ y/t. 

In conclusion, we have analyzed the problem of grain 
growth in a vortex polycrystal studying the dynamics 
of grain boundaries. We have obtained estimates for 
the grain size that compare well with experiments on 
NbMo crystals, and derived the law for thermally acti- 
vated grain growth. In general, our theory should apply 
to field cooling experiments in which there is a competi- 
tion between the ordering stress and the pinning stress. 
The resulting polycrystalline structure, at least in the 
weak pinning regime, represents a metastable state in 
which the system is trapped during its evolution towards 
the stable Bragg glass phase. The possibility of a thermo- 
dynamically stable polycrystalline vortex state, recently 
suggested by experiments 0,03, s * m remains to be con- 
firmed theoretically. 
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